Reentrant charge order transition in the extended Hubbard model 



R. Pietig, R. Bulla and S. Blawid 

Max- Planck- Instttut fur Physik komplexer Systeme, Nothmtzer Str. 38, 01187 Dresden, Germany 

We study the extended Hubbard model with both on-site and nearest neighbor Coulomb repulsion 
(U and V, respectively) in the Dynamical Mean Field theory. At quarter filling, the model shows 
a transition to a charge ordered phase with different sublattice occupancies nA^riB- The effective 
mass increases drastically at the critical V and a pseudo-gap opens in the single-particle spectral 
function for higher values of V. The V c (T)-curve has a negative slope for small temperatures, i.e. 
the charge ordering transition can be driven by increasing the temperature. This is due to the higher 
spin-entropy of the charge ordered phase. 
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The possibility of crystallization of electrons due to 
their long-range Coulomb repulsion was first proposed 
by Wigner 0. He considered an electron system in a 
uniform positive background at sufficiently low densities. 
The Wigner lattice is formed when the gain in Coulomb 
energy due to the localization of the electrons exceeds the 
gain in kinetic energy for the homogeneous electron dis- 
tribution, ft is experimentally realized in the two dimen- 
sional electron gas in a GaAs/AlGaAs heterostructure 
g . Due to the reduced dimensionality, the effect of the 
Coulomb interaction is enhanced so that the transition 
to the ordered state occurs at experimentally accessible 
electron densities. 

Crystallization of charge carriers (charge ordering) can 
also be observed in three dimensional systems, even at 
very high densities Here the kinetic energy of the 
electrons or holes has to be reduced drastically for the 
charge ordered state to become possible. In 4f-electron 
systems it is the small hybridization of the well localized 
4f-orbitals that leads to a reduced kinetic energy. An 
example is Yb4As3 where a first order charge ordering 
transition occurs at T c w 295K J|||. The carrier con- 
centration in Yb4As3 (approximately one hole per four 
Yb ions) is considerably larger than typical values for a 
Wigner lattice. The kinetic energy of the electrons can 
also be reduced by the interaction with lattice and spin 
degrees of freedom. An interplay of these mechanisms is 
responsible for the charge order transition occurring in a 
variety of rare earth manganites (e.g., in Lai_ x Ca x Mn03 
for x>0.5 §). 

In all examples mentioned so far, the charge ordered 
phase is the ground state. However, a melting of 
the charge ordered state on decreasing the temperature 
(i.e. a reentrant transition) has been found recently in 
Pro.65(Cao.7Sro.3)o.35Mn0 3 § and in LaSr 2 Mn 2 7 g. 

In this Letter, we investigate the simplest model which 
allows for a charge ordering transition due to the com- 
petition between kinetic and Coulomb energy. The ex- 
tended Hubbard model M 

H = 1 E ( C L C JV + C ]a C i°) ™ ^ E 4a C ^ 
<ij>a ia 



i <ij> 

describes fermions on a lattice with an on-site Coulomb 
repulsion U, a nearest neighbor Coulomb repulsion V 
and a hopping matrix element t. The c\ a (ci a ) denote 
creation (annihilation) operators for a fermion at site i 
with spin a, the rij are defined as ni = + where 
riio- = c\ a Ci a and Yl<ij> indicates the sum over nearest 
neighbors. 

In the following we study the extended Hubbard 
model eq. (f|) within the Dynamical Mean Field Theory 
(DMFT) 00, i- e - in the limit of infinite lattice coordi- 
nation number z. In order to define a nontrivial limit as 
z^oo, the parameters t and V are rescaled as t — > t/\fz 
and V — > 2V/z. This leads to a drastic simplification of 
the self-energy diagrams. The self-energy becomes local 
and in particular the nonlocal Coulomb term contributes 
only at Hartree level, i.e. the V-teim simply acts as a shift 
of the chemical potential E|. Therefore, the hamiltonian 

ia j 

+ 1 E (4r c J<r + 4^) +U Y1 4l c ^4l c ii> ( 2 ) 
<ij>a i 

(Y^'j indicates the sum over the nearest neighbors of i). 
leads to the same z — > oo limit as eq. ([!]) after rescaling t 
and V. 

As we are interested in charge ordered phases with dif- 
ferent occupancies on the two sublattices A and B of a 
bipartite lattice, we have to generalize the DMFT equa- 
tions to allow for solutions with ^ fiB ■ For this, the 
model (||) is mapped self-consistently on two Anderson 
impurity models (one for sublattice A and one for sub- 
lattice B). In the Bcthe lattice case, where the density of 
states is given by D(e) — -^n-\/W 2 — e 2 (we set W = l as 
the unit for the energy scale), the self-consistency equa- 
tions simplify to the form: 

e% + A A (iu) n ) = 2Vn B - Li+~G B (ioJ n ), (3) 
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ef + A B (««/») = 2Vn A - /i + -G A (^„), (4) 

(Ga/b(^ r ) are the Green functions for the A/B sublat- 
tice and Aa/b^^™) are the hybridization functions be- 
tween impurity and effective conduction band; Ua/b de- 
note the sublattice occupancies and are the on-site 
energies of the effective impurities). 

The remaining problem is the iterative solution (i.e. 
the calculation of the Greens functions) of an effective 
single impurity Anderson model. We use three differ- 
ent methods, an Exact Diagonalization technique (ED) 
for finite temperatures, the Non-Crossing Approxima- 
tion (NCA) and the Numerical Renormalization Group 
method (NRG). The ED method diagonalizes an impu- 
rity model with a finite number N of conduction band 
orbitals. Within the self-consistency procedure we de- 
fine the mapping of the full Green functions Ga/b(*w„) 
to the hybridization functions A b /a(* w k) by expanding 
both sides of @,(Q) in powers of (iLUn)^ 1 and match co- 
efficients up to order (iuj n )~ 2N . This approximation is 
similar in spirit to the projection method based on the 
continued fraction representation used for T = |ll]] . The 
NRG is applied here for the first time to a particle-hole 
asymmetric problem within the DMFT. The method is 
an extension of earlier work on the (not extended) sym- 
metric Hubbard model HjJ. Details of the NCA approach 
are summarized in 

Fig. shows the (T-V)-phase diagram for U — 2 and 
quarter filling. The results of the different methods agree 
remarkably well in their corresponding range of applica- 
bility. For high temperatures (T>0.4), the ED results 
for N = 5 shown in Fig. [l] can already be obtained from an 
N = 3 calculation (within numerical accuracy). The ED 
method cannot be used for very low temperatures since 
only a small number of orbitals is taken into account. The 
NCA is applicable down to much lower temperatures and 
we find that the slope of the V c (T)-curve changes its sign 
at T « 0.1. The NCA encounters problems in the very 
low temperature limit. Nevertheless, the extrapolation 
of the V c (T)-curve to T = agrees well with the critical 
V obtained from the NRG method V C (T = 0)«0.66. 

The reentrant behavior in the interval 0.47< V<0.66 
is due to the higher spin entropy of the charge ordered 
state. As we have neglected the possibility of an addi- 
tional magnetic ordering in the charge ordered state, the 
contribution of the spin degrees of freedom to the entropy 
is approximately ln(2) per occupied site. The entropy of 
the homogeneous phase, which is a Fermi liquid, grows 
linearly with T. On further increasing the temperature, 
the increasing charge entropy of the homogeneous phase 
dominates and the system shows the usual melting be- 
havior. 

The inset of Fig. [l] shows the ED-result for the In- 
dependence of the sublattice occupancies tt-a and ub 
(T = 0.2, U = 2). The transition is clearly continuous, 
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FIG. 1. Phase diagram for the extended Hubbard model 
(U — 2, quarter filling). The various symbols show the results 
for the phase boundary between the homogeneous and the 
charge ordered phase (CO) from ED (circles), NCA (squares) 
and NRG calculations (diamond). The inset shows the ED re- 
sult for the Independence of the lattice occupancies for (7 = 2 
and T = 0.2. 

in contrast to the result for T = where the NRG gives a 
first order phase transition with a jump in the order pa- 
rameter 5n — nA — n-B from to Sn^ 0.7. Unfortunately, 
it is not possible to clarify numerically how the first order 
transition at T = evolves from the continuous transition 
at finite T . There are numerical indications that 5n in- 
creases more rapidly at the transition when T is reduced. 
However, the convergence of the iterative procedure is ex- 
tremely slow near the phase boundary. 

The NRG results for the A and B spectral functions 
for T = are shown in Fig. ^. Below V = V c , the A 
and B spectral functions are equal and independent of V 
(the extended Hubbard model in the DMFT reduces to 
the ordinary Hubbard model as long as the homogeneous 
phase is considered; the solution shown in Fig. § for V — 
0.65 is therefore also the solution for the Hubbard model 
with V = 0) . The discontinuous transition to a finite 5n 
at V c is reflected in a redistribution of spectral weight in 
the spectral functions above V c . 

For V = 0.665 the system is still a Fermi liquid with 
a strongly enhanced effective mass (m*/m ss 120) and a 
very narrow quasiparticle resonance at the Fermi level. 
The effective mass increases further on increasing V and 
an extreme narrow quasiparticle resonance is still present 
for V= 1.5 (the width of this resonance is below the en- 
ergy resolution used in Fig. so that the peak is not 
visible). Although the system is still a Fermi liquid on 
an extremely small energy scale, a pseudo-gap in the A 
and B spectral functions gradually develops at V~0.9. 
The spectral weight within the pseudo-gap decreases ap- 
proximately as 1/V 2 for large V. 

The NCA-result for the spectral function for U = 2, 
T = 0.2 and various values of V is shown in Fig. ||. The 
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FIG. 2. NRG-results for the A and B spectral functions 
p(u) (solid and dashed lines, respectively) for U — 2, T — 
and different values of V. Below the critical V c ~ 0.66, the 
spectral functions of both sublattices are equal. 
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FIG. 3. NCA-results for the A and B spectral functions 
p(aj) (solid and dashed lines, respectively) (7 = 2, T = 0.2 and 
different values of V. For V>0.58 charge order sets in and a 
pseudo-gap develops at the Fermi level. 



very narrow quasiparticle peak seen at T = is absent 
because the temperature exceeds the low energy scale 
associated with the Fermi liquid behavior. In contrast 
to the T = 0-case, the transition is continuous which is 
reflected in the continuous transfer of spectral weight at 
the transition. A comparison of the NRG and NCA re- 
sults for V — 0.9 and V — 1.5 shows a good agreement 
of the main features of the spectral functions. However, 
the resolution of the high energy features is limited in 
the NRG since the spectra are obtained by broadening a 
discrete set of <5-peaks. 

In the homogeneous phase, the NCA spectral function 
shows the splitting in upper and lower Hubbard band. 
The weight of the upper Hubbard band in the B-spectral 
function is considerably larger compared to the NRG re- 
sult for T=0. The weight decreases on increasing V in 
the charge ordered phase because the occupancy of the B- 
sites is suppressed. In addition, the bands in the spectral 
function narrow due to the reduced hopping of the elec- 
trons in the charge ordered phase. In the limit of large 
V, perfect charge order evolves (jia — > 1 and ne — > 0). 
Therefore, in this limit the total spectral function of the 
extended Hubbard model has a simple three-peak struc- 
ture with peaks at E\ = —ji, E% = U — /i and £3 = 2V — \i. 
The peak at E\ corresponds to the process of taking out 
one electron from a (singly occupied) A-site, the peak at 
E2 is due to the addition of an electron to an A-site and 
the peak at E3 comes from the addition of an electron to 
a B-site. 

It is tempting to explain the opening of the charge- 
order gap by the change of the unit cell alone, reflecting 



the broken AB-symmetry of the lattice. Although the 
AB-symmetry is of course broken at the charge order 
transition, our results for the spectral function cannot 
be understood within a one-particle picture. First of all, 
the charge-order gap visible in the spectral function is 
not a real hard gap but only a pseudo-gap (i.e. an energy- 
regime of strongly reduced spectral weight) even at zero 
temperature. From an experimental point of view this 
difference is irrelevant. However, a pseudo-gap cannot 
be obtained by a simple band splitting due to a reduced 
symmetry. To demonstrate the many-particle charac- 
ter of the bands, we have doped the extended Hubbard 
model away from quarter-filling (see Fig. Note that 
within our approach the charge ordered state on the A/B- 
sublattice exists also for n^l. A rapid weight transfer- 
between the high- and the low-energy scale is observed, 
which is a characteristic feature of strongly correlated 
electron systems. A well examined example is the usual 
half filled Hubbard model, where the weight of the lower 
Hubbard band decreases faster than the doping x ]l5],[l6| . 
In the extended Hubbard model we find the same behav- 
ior. 

In conclusion, we have found that the charge order 
transition in the infinite dimensional extended Hubbard 
model shows a number of very interesting features. For 
low temperatures, the system is a Fermi liquid for both 
V < V c and V > V c . Nevertheless, the effective mass 
undergoes a rapid change at the transition. The corre- 
sponding very narrow quasiparticle peak at the Fermi 
energy will hardly be seen in photoemission experiments 
of charge ordered materials. However, the huge mass en- 
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FIG. 4. NRG-results for the A and B spectral functions 
p(uj) (solid and dashed lines, respectively) for V — 1.0, T = 
and different occupancies n. A transfer of weight from the 
lower Hubbard band (of both A and B spectral functions) to 
the Fermi level is clearly visible. 



hancement at the transition has drastic consequences on 
the transport properties, e.g. on the resistivity. This is 
the case in YD4AS3, where the charge order transition is 
accompanied by a jump in the resistivity. 

In our model calculations, we find a real first order 
transition only for T = 0, although the transition remains 
sharp also for small T^O. The occurrence of a first order 
transition at finite temperatures in the above mentioned 
examples indicates that lattice degrees of freedom are 
involved. Indeed, the charge order transition in Yb4As 3 
is accompanied by a structural phase transition. 

For a particular range of the nearest neighbor Coulomb 
repulsion V the extended Hubbard model shows reen- 
trant behavior, i.e. a transition from the charge or- 
dered to the homogeneous phase with decreasing tem- 
perature. This behavior is experimentally seen in 
Pro.65(Cao.7Sro.3)o.35Mn0 3 § and in LaSr 2 Mn 2 7 @. 
Moreover, the bandwidth in Pro.65(Cai_ y Sr y )o.35Mn03 
can be varied by doping with Sr or by applying a mag- 
netic field. For sufficiently small kinetic energy of the 
electrons only one charge order transition occurs. In- 
creasing the kinetic energy, a reentrant regime is found 
until for high dopings or magnetic fields the charge order 
is totally suppressed. Although the localization mech- 
anisms in the rare-earth manganites are more complex 
than the one induced by a nearest neighbor Coulomb 
repulsion, we find the same behavior in our model as 
function of V/W. We take this as an indication that the 
occurrence of a reentrant charge order transition does not 
depend on details of the localization mechanism. 

Finally, in photoemission experiments of charge or- 



dered materials the many-particle character of the sys- 
tem should be seen in a weight-transfer upon doping with 
additional charges. 

The authors would like to thank P. van Dongen and P. 
Thalmeier for helpful discussions. 
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